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It is pointed out that the "exact" solutions given in the recent Comment by A. Rabhi on the 
paper [N. Dinh Dang "Energies of the ground state and first excited 0"'' state in an exactly solvable 
pairing model", Eur. Phys. Jour. A 16 (2003) 181] do not correspond to the standard exact 
solution of the well-known two-level pairing model. Other issues raised in the Comment reiterate 
the discussions already published in the original paper by N. Dinh Dang. 
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In Ref. several approximations, namely the BCS approximation, Lipkin - Nogami method, random-phase ap- 
proximation (RPA), quasiparticle RPA (QRPA), the renormahzed RPA (RRPA), and renormalized QRPA (RQRPA), 
are tested by calculating the ground-state energy and the energy of the first excited 0+ state using the well-known 



I exactly solvable model with two ponmetric levels interacting via a pairing force. The author of the recent Comment 

. |5| discusses three points of Ref. namely (A) the boson mapping in Sec. 4.1.1, (B) the approximations leading to 

^ ' the RPA matrices (58) - (61) and those obtained in Ref. and (C) the approximation (75). 

^ ■ (A) The author of 0] claimed that the one-boson energy w^py^ given by Eq. (53) of [3 is erroneous. As the proof he 

OO ■ introduces two "exact" solutions for the 0+ energy in Fig. 1 of Q]- However, these "exact" solutions are completely 

• different from the standard exact solution of the well-known two-level model under consideration, which has been 

^ ■ obtained using the SU(2) algebra in many papers \See Sec. 2 and Fig. 3 of Ref. 0, or Refs. 0,01, e.g.]. In particular, 

, as has been mentioned at the end of Sec. 4.1.1 of the solution given by Eq. (53) of 0] is exactly the same as that 

J — ■ obtained for the first time by Hogaasen-Fledman in Ref. 3 using the space-variable technique. The author of the 

f — , Comment 0, however, fails to reproduce the boson and exact solutions by Hogaasen-Feldman, saying that he does 

■ not understand it ^. 

, (B) Sec. 4.1.2 just discusses two approximations (62) and (63). The former is based on the exact commutation 

■ relations in Eq. (10) and leads to the matrices (58) - (61). The latter leads to those in Ref. The difference is the 
^0 ' factor of 2 in the denominator of the second term at the right-hand side of Eq. (58) in [ij and the corresponding term 
■>«.^ in Eq. (38) of The omission of q-teim was considered in [l| also as a possible approximation, which leads to the 

' appearance of the spurious mode as has been pointed out in 0] and repeated by the author of It is worthwhile 

I , to study this approximation since there have been numerical calculations within the RPA neglecting the so-called 

^ ' scattering terms, which have the same origin as that of the q-term considered here (See the discussion in (b) on page 

^ 185 of [l|). In such calculations the parameters of the effective interaction are usually readjusted in such a way that 

^ ■ that the energy of the spurious mode is zero to compensate for such effect [See, e.g. Ref. 0]]. To my knowledge, 

^ I there is no exact way to take into account the scattering term within the QRPA so far. Approximations to take into 

• ^ . account the scattering term lead to the extended QRPA 7*] or modified QRPA ^ . 

' (C) The approximation (75) in ,1] has been introduced so that one can compare the exact solution, the phonon 

^ , solution obtained within the fermion formalism in Sec. 4.1.2 with that obtained within the one-boson mapping in 

■ - - ' Sec. 4.1.1. The conclusion is that the exact solution can be approximately considered as a mixture of superfluid and 
normals fluid states as developed in Sec. 4.3. In this sense the approximation (75) mixes iV ± 2 states. 
The RPA operator for the additional mode in the present two-level model is 

Qi = - Yi4 ■ (1) 

The one for the removal mode is 

Ql = XiA\ - Y2AI . (2) 

It is clear that neither Q\ nor Q\ conserves the particle number on each level because, according to the exact 
commutation relations in Eq. (5) of I*], they do not commute with Ni (i =1,2), except for Xi = O 01 Yi — 0. So 
doesn't the operator Q'^ = Q\ + Q\ in Eq. (75) except for Xi = 0. However, in the boson formalism, based on 



^ The derivation in Eqs. (2) - (10) in the Comment Q was actually sent by Dinh Dang to Rabhi after the latter failed to reproduce the 
result of Eq. (53) of Qj as well as that obtained by Hogaasen in |^ 
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the mapping (49) and (50) of Q with bi = b2 = b as has been discussed in the Comment [3, one obtains the boson 
images of these commutation relations as follows 

[Ni,Qi]b = [2n - 2btb, X2h^ - Yib] = -2X2bt - 2Yib , (3a) 

[N2,Q% - [2btb, X2bt - Yib] = 2X2bt + 2rib , (3b) 

[Ni,Q% - [2r! - 2btb, Xib - Fabt] = 2Xih + 2Y2h^ , (3c) 

[7V2, Ql]b = [2btb, Xib - Fab^] = -2Xib - 2y2b^ , (3d) 
Summing up Eqs. (|3a|l and Ij3bp . one finds 

[iVi+iV2,Ql]b = [A^,Ql]b = 0. (4) 

Summing up Eqs. (|5c|l and (|3d|l . one finds 

[Ni+N2,Ql]b^[N,Ql]i,^0 . (5) 

Hence 

[iv,Qt]^[iV,gt + g,] = o. (6) 

These results show that, in the boson formalism, although the additional and removal operators do not conserve the 
particle number Ni on each level, they do conserve the total particle number N = Ni + N2. So does the operator 
gt ^ Qt + Qt because of Eq. 

In conclusion, the present Reply clarified several issues raised in the recent Comment 0] on the paper Q regarding 
the RPA and QRPA for an exactly solvable pairing model. 
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